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Abstract 

Human  Brains  are  highly  convoluted  surfaces  with  multiple  folds.  To  characterize  the 
complexity  of  these  folds  and  their  relationship  with  neurological  and  psychiatric  conditions, 
different  techniques  have  been  developed  to  quantify  the  folding  patterns  and  gyrification  of  the 
brain.  In  this  work,  we  propose  a  new  geometric  approach  to  measure  the  local  gyrification  of 
human  brains  from  magnetic  resonance  images  (MRI).  This  approach  is  based  on  intrinsic  3D 
measurements  that  relate  the  local  brain  surface  area  to  the  corresponding  area  of  a  tightly 
wrapped  sheet.  These  quantities  are  efficiently  and  accurately  computed  solving  geometric 
partial  differential  equations.  The  presentation  of  the  geometric  framework  is  complemented 
with  experimental  results  for  brain  complexity  in  typically  developing  children  and  adolescents. 
Using  this  novel  approach,  we  provide  evidence  for  developmental  alterations  in  brain  surface 
complexity  throughout  childhood  and  adolescence. 


Keywords:  Gyrification  index,  cortical  complexity,  gyri,  sulci,  age  population  studies,  brain 


evolution 
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1  lotroductian 

The  conical  surface  of  a  human  brain  evolves  from  a  smooth,  or  llsencephalic  surface,  to 
a  highly  convoluted  surface  during  the  third  irimesierof  feial  life  (Welker.  1990).  This  process  is 
known  as  gyrificaiion,  and  by  the  lime  of  binh  the  human  brain,  although  smaller,  has  an 
appearance  very  much  like  an  aduli  brain.  Even  though  a  conical  surface  is  considered  to  be  a 
topological  sphere,  the  gyri  and  sulci  that  form  the  fissures  and  folds  can  be  complicated  (Welker. 
1990).  It  has  been  suggested  that  a  higher  degree  of  folding  of  the  conical  surface  indicates  a 
progressive  evolution  of  cortical  complexity  in  humans  (Zille.s  et  al.  1988).  The  increased 
conical  folding  associated  with  human  phylogeny  has  produced  a  highly  efficient  mechanism  to 
increase  cortical  gray  matter  while  optimizing  a  smaller  overall  brain  size.  Cortical  complexity, 
which  is  considered  a  measure  of  gyrification.  reveals  the  general  organization  in  the  human 
cortex  and  the  inherent  structural  configuration  of  the  brain. 

It  is  of  great  interest  to  characterize  the  general  brain  shape  and  provide  valu^le 
information  associated  with  pathology  through  quantitative  measures.  The  cortical  complexity 
has  been  used  to  evaluate  abnormalities  of  the  brain's  surface  morphology  in  various 
neurological  and  psychiatric  conditions  and  in  disorders  of  cognition  (White  and  Hilgetag.  2(X)8). 
Estimates  of  coitical  complexity  derived  from  the  number  of  voxels  which  occupy  the  boundary 
between  the  gray  matter  and  CSF  in  the  prefrontal  cortex  have  shown  that  abnormalities,  such  as 
smaller  prefrontal  volume,  greater  sulcal  width  and  gyral  pattern  abnormalities,  are  present  in 
first-episode  schizophrenia  (Wiegand  et  al.  200S).  Since  there  is  a  proposed  connection  between 
development  and  gyrification  (Armstrong  et  al.  1995),  and  between  gyrification  and  connectivity 
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(Van  Essen.  19^).  these  findings  may  help  ascertain  the  underlying  neurobiology  of 
schizophrenia  and  other  psychiatric  and  neurological  disorders. 

There  is  growing  interest  in  how  the  human  brain  develops  gyri  and  sulci>  how  the 
differences  in  the  degree  of  convolutions  affect  brain  function,  and  how  the  cortical  folding 
evolves  and  whether  this  evolution  is  related  to  human  age>  gender,  ethnic  group  or  health  status. 
These  mysteries  of  human  brain  development  have  prompted  researchers  to  develop 
measurements  to  quantify  the  degree  of  complexity  of  the  cortical  surface.  The  gyrification  index 
(GI)  proposed  by  Zilles  et  ai.  (Zilles  et  al.  19S8)  was  developed  to  quantify  the  degree  of  cortical 
folding.  Their  two  dimensional  global  measurement  used  an  estimate  of  the  ratio  between  the 
total  outline  (including  sulcal  folds)  and  the  superficially  exposed  outline  using  coronal  sections 
of  postmortem  brains.  This  mea.sure  ha.s  been  applied  to  study  both  the  phylogeny  (Zilles  et  al. 
1988)  and  ontogeny  (Armstrong  et  al.  1995)  of  cortical  gyrification.  Brains  that  have  higher 
degrees  of  conical  folding  (i.e..  increased  cortical  complexity)  yield  larger  values  of  this  GI. 

Anterior  to  posterior  maps  of  human  GI  have  shown  greater  gyrification  in  the  frontal, 
temporal,  and  parietal  lobes  of  the  brain  (Zilles  et  al.  I9SS).  However,  due  to  the  oblique 
orientation  and  manual  delineation  (Schaeret  al.  2(X)8;  Zilles  et  al.  1988).  the  intrinsic  3D  nature 
of  the  brain  surface  is  not  taken  into  account  and  this  gyrification  index  does  not  completely 
avoid  a  biased  estimate.  For  example,  the  GI  may  be  altered  if  (he  slice  orientation  is  slightly 
different  within  the  same  subject.  Since  surface  morphology  can  be  considered  a  two> 
dimensional  (2D)  surface  in  three-dimensional  (3D)  space,  it  is  important  to  consider  a  full  3D 
gyrification  index  to  eliminate  the  shortcomings  of  the  coronal  2D  approach.  A  number  of 
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neuroimagiRg  software  packages  are  currenUy  available  lo  generate  3D  reconsiruciions  of  human 
brains  (eg.  Freesurfer  (Dale  el  al.  1999) SurfRelax  (Larsson,  200!)^.  and  BrainVisa  (Han  ei  al, 
2004)^),  These  software  packages  provide  the  early  preprocessing  stages  necessary  to  develop  a 
3D  gyrification  Index,  and  one  such  3D  technique  has  been  recently  developed  (Schaer  etal. 
2008),  following  in  part  computational  techniques  developed  in  (Kao  ei  al,,  2(XI7), 

Previous  studies  have  proposed  reliable  computational  algorithms  to  quantify  gyrification 
asing  either  fully  automatic  or  semi-automatic  approaches  (Magnotta  et  al.  1999;  Schaer  et  al, 
2008;  Zilles  et  al.  1988).  The  manually  generated  landmark  studies  mentioned  above  defined  GI 
as  the  ratio  of  lengths  in  2D  coronal  slices  (Zilles  et  al,  1988)  and  produced  a  GI  value  for  each 
coronal  slice  of  each  hemisphere.  A  curvature-based  approach  introduced  in  (Luders  et  al.  2{XI6; 
Schaer  et  al.  2008)  evaluates  the  smoothed  absolute  mean  curvature  on  parameterized  conical 
mesh  models,  revealing  that  females  have  increased  gyrification  in  frontal,  parietal,  temporal, 
and  occipital  conices  compared  to  males.  This  curvature  based  approach  found  differences 
which  were  not  previously  identified  using  the  traditional  GI  approach.  More  recently,  a 
localized  3D  gyrification  index  (IGI)  (Schaer  et  al.  2008)  was  developed  and  applied  to  a  group 
of  children  affected  by  22ql  1,2  Deletion  Syndrome.  This  technique  used  3D  triangulated  mesh 
reconstructions  of  cortical  surfaces  and  outer  hull  surfaces,  and  measured  the  amount  of  cortical 
surface  buried  in  the  sulci  by  constructing  and  relating  this  to  a  non-intrinsic  sphere  with 
different  radii.  The  authors  used  the  approach  developed  by  Kao  et  al,  (Kao  et  al,  2007)  to 

’  iee  hllp  //surFcr.nmr.tngn.harvard.cdu/ 

‘  Suiffitkix,  sac  htlo://www.cns.nyij.frdij/~|onas/sQflwarc.lilml 
'  SreinVIso,  see  lillp://brainvlsa  Info/ 
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generate  ibe  outer  hull  surface  and  the  id  for  each  point  on  the  conical  surface  Is  obtained 
through  a  depth-weighted  sum  of  neighboring  points.  We  wiJJ  comment  on  this  approach  after 
we  introduce  our  fully  geometric  G1  measure. 

In  this  paper,  we  build  upon  our  previous  work  (Kao  et  al.  2007)  and  propose  a  3D 
geometric  approach  for  the  automatic  compulation  of  global  and  local  gyrification  indices  of  the 
human  brain.  By  finding  a  geometrically  corre.sponding  region  on  the  outer  hull  surface  for  any 
selected  region  of  interest  on  the  conical  surface,  we  define  a  3D  local  GI  as  the  area  ratio 
between  the  selected  region  and  the  corresponding  region  of  a  tightly  wrapped  sheet. 
Furthermore,  our  3D  Gl  can  be  weighted  by  local  quantities,  i.e.,  curvature  and  sulcal  depth,  and 
is  fully  intrinsic  and  different  than  the  method  proposed  by  Schaer  et  al.  (Schaer  et  al,  2008)  in 
that  it  does  not  depend  on  a  chosen  radius  and  a  corresponding  non  intrinsic  ball  to  determine  the 
region  used  to  calculate  the  IGI.  The  incorporation  of  the  robust  sulcal  depth  computation 
developed  in  (Kao  et  al.  2007)  as  part  of  the  GI  measurement  is  important  to  characterize 
different  levels  of  convolutions  in  human  brains. 

In  Section  2.  we  introduce  the  proposed  novel  algorithms  to  calculate  3D  gyrification 
indexes.  Quantitative  results  for  our  proposed  gyrification  indexes,  and  comparison  with 
previous  approaches,  are  presented  in  Section  3.  For  this,  we  present  the  application  of  our 
algorithm  to  a  population  of  typically  developing  children.  Finally,  we  conclude  the  paper  in 
Section  4  briefly  de.scribing  our  contributions  and  findings  and  proposing  further  lines  of  study  in 
the  area  of  computational  approaches  for  measuring  GI, 
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2  MateriaJs  and  Methods 


A  flow  chart  explaining  the  main  steps  of  the  proposed  2D  geometric  GI  computation  algorithm 
is  shown  in  Figure  I,  First,  a  T1 -weighted  MRI  volume  is  segmented  into  different  regions,  e  g., 
white  matter  (WM),  gray  matter  (GM).  and  cerebral-spinal  fluid  (CSF)  regions,  u.sing  for 
example  freesurfer  (Dale  el  al.  1999).  We  consider  the  cortical  surface  represented  as  a  triangular 
mesh  as  the  input  to  our  algorithm.  Second,  we  construct  an  outer  hull  surface  covering  the  sulcal 
regions  in  a  shrink-wrap  fashion.  Third,  we  define  the  geodesic  depth  a.s  the  shortest  distance 
connecting  each  point  on  the  cortical  pial  surface  to  a  point  on  the  outer  hull  surface,  the 
connecting  path  staying  inside  the  sulcal  regions,  and  compute  the  depth  by  using  an  efficient 
fast  sweeping  algorithm.  Fourth,  a  patch  on  the  cortical  surface  is  selected  and  its  boundary  is 
extracted.  Fifth,  the  corresponding  part  of  the  selected  patch  on  the  outer  hull  surface  is 
computed.  Sixth  and  finally,  we  define  the  local  3D  GI  as  the  ratio  between  the  area  of  the  patch 
on  the  cortical  surface  and  the  area  of  the  corresponding  patch  on  the  outer  hull  surface.  This 
ratio  can  be  weighted  by  the  local  sulcal  depth.  We  now  explain  each  one  of  these  steps  in  detail. 


2.1  Surface  Extraction  and  Depth  Computation 


The  input  of  our  algorithm  is  a  T I -weighted  MRI  human  brain  volume.  For  skull 


stripping,  segmentation  of  the  brain  into  white  matter,  gray  matter,  and  cerebrospinal  fluids,  and 
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exiraciion  of  a  topologically  correct  triangular  me^h  representing  the  piaJ  (GM*CSF)  surface  of 
the  cerebral  corieXs  we  use  in  (his  work  the  publicly  available  software  freesurfer  (Dale  ei  al, 

1999).  By  we  denote  a  triangle  mesh  with  faces  T, . T, ;  we  require  to  be  a  closed  and 

orientable  2-manifold  in  Euclidean  2*space,  Then  we  use  a  regular  grid  to  derive  an  implicit 
representation  of  the  pial  surface  by  computing  the  signed  distance  function  to  the  surface  on  a 
Cartesian  grid  (Kao  et  al.  2007).  By  using  a  level  set  technique  we  compute  an  outer  hull  surface 
that  encloses  the  pial  surface  in  a  shrink-wrap  type  fashion.  The  outer  hull  i.s  such  that  one  can 
still  distinguish  the  gyri>  but  the  sulcal  regions  are  now  “filled."  Then  we  define  a  depth  measure 
for  the  pial  surface  as  (be  shortest  distance  that  connects  each  surface  point  to  a  point  on  the 
outer  hull  surface  such  that  the  connecting  path  remains  inside  the  sulcal  regions.  The 
computational  realization  uses  an  efficient  fast  sweeping  algorithm  (Kao  et  al.  2005;  Tsai  et  al. 
2002;  Zhao.  2005). 

As  shown  in  Figure  1(a),  once  we  have  the  explicit  form  of  the  pial  surface,  we 
compute  the  signed  distance  function  to  the  surface  on  a  Cartesian  grid.  In  the  implicit  form,  the 
pial  surface  becomes  the  zero  level  set  (cp  »  0}  of  the  signed  distance  function  r'  -*  R 
(Osherand  Sethian.  1988).  The  details  for  retaining  this  signed  distance  function  can  be  found  in 
(Kao  et  al,  2007). 

After  retaining  both  an  explicit  (riangular-mesh  representation  and  an  implicit  level-set 
representation  on  a  Cartesian  grid,  we  compute  the  outer  hull  surface  using  a  morphological 
closing  operation  applied  to  the  level  set  function  ®  (Osher  and  Sethian.  1988).  That  is,  we  first 
move  the  pial  surface  outward  by  a  time  parameter  T  .  and  then  we  move  the  surface  inward  by 
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ihe  same  amount  of  lime.  In  our  algorithm,  we  choose  r  •  lo  (mm/unii  time),  since  T  needs  to 
be  large  enough  to  fil I/close  (he  sulcal  regions  and  also  small  enough  to  keep  the  overall  shape  of 
the  gyri.  The  governing  equation  is 

®,  ♦  V(/)J  V®  la  0 
®  l>.  v.j.O)  ■  ®  U,  y.r), 


where 


for  /  <  r 
for  T<i  <1T 


Thus,  (he  impUci(  representation  of  the  outer  hull  surface  is  given  by 


‘♦'(.a)  smmj  ®  ( t.  y.2,2T),®(A.  v,:.0)|. 


This  minimization  guarantees  that  (he  ou(erhull  surface  covers,  bu(  doe.s  no(  pene(ra(e  (he  pial 
surface.  In  Figure  1(b)  we  show  (he  compu(ed  ou(er  hull  surface  for  (he  given  piaJ  surface  shown 
in  Figure  1(a). 

Once  we  generare  (he  ou(er  hull  surface,  we  apply  (he  fast  sweeping  method  (Kaoet  al. 
2005;  Tsai  et  al,  2003:  Zhao.  2005)  to  calculate  the  geodesic  depth  for  points  on  the  pial  surface. 
The  geodesic  depth  we  define  corresponds  to  the  shortest  path  from  (he  given  pial  surface  point 
to  (he  computed  outer  hull  which  do  not  cross  the  surface  of  the  brain  (Kao  et  al.  2007).  and  thus 
we  apply  our  distance  computation  algorithm  to  the  restricted  (CSF)  region  between  Ihe  outer 
hull  and  (he  pial  surface  I'l'^Oand  ®a0}.ln  Figure  )(c).  we  show  a  top  view  of  (he  bottom 
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pan  of  the  piaJ  surface.  The  color  coding  corresponds  to  ibe  computed  geodesic  depth  where 
shallow  regions  are  colored  in  blue  and  deep  regions  are  colored  in  red. 


2.2  ('orre.sponding  Region  ad  (he  Outer  Hull  Surface  ComputatiAn 


In  this  stage,  we  first  select  a  region  of  intere.si  (ROI)  on  the  pial  surface.  The  region  can 
be  chosen  using  several  different  algorithms  and  can  involve  any  sulcus  or  gyrus  of  interest. 
Alternatively,  the  region  can  be  restricted  by  implementing  either  a  depth  or  curvature  threshold. 
After  obtaining  the  ROI  on  the  pial  surface,  the  corresponding  region  on  the  hull  surface  is 
computed  in  two  steps. 

First,  the  boundary  of  the  pial  surface  for  the  selected  region  is  extracted.  Then  we 
sample  points  on  this  boundary  and  find  the  boundary  points  for  the  corresponding  region  on  the 
computed  outer  hull  surface.  This  is  done  following  the  negative  gradient  of  the  previously 
computed  sulcal  depth.  We  extract  the  path  s  for  each  starting  point  situated  on  the  boundary  of 
the  selected  pial  region  by  solving  tbe  following  equations,  and  stop  when  we  reach  the  outer 


dY 

-  ■  -Vrf 

X(0)  S 
X{t)  m  Bl 


hull  surface: 
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where  d  is  ibe  geode?>ic  depth  we  computed  in  Section  lA,  U  a  point  on  the  l>oiindary  of  the 
selected  region  on  the  pial  surface,  and  5.'  is  then  the  corresponding  boundary  point  on  the  outer 
hull  surface.  In  this  way  we  obtain  all  the  corresponding  boundary  points  on  the  outer  hull 
surface,  in  an  intrinsic  and  geometric  fa.shion,  and  need  to  close  the  boundary.  We  connect  the 
computed  outer  hull  boundary  points  by  using  the  shortest  distance  along  the  edges  of  the 
triangles.  An  alternative  choice  would  be  to  use  the  geodesic  distance  on  the  triangular  mesh 
(Peyre  and  Cohen,  2006),  Since  is  a  closed  and  orientable  2- manifold  in  Euclidean  3-space, 
we  can  march  inside  it  and  find  the  triangles  restricted  in  the  constructed  boundary  on  the  outer 
hull  surface.  The  region  defined  by  these  restricted  triangles  corresponds,  on  the  outer  hull 
surface,  to  the  selected  pial  region  of  interest.  In  Figure  1(d),  we  show  a  selected  region  in  red  on 
the  pial  surface,  and  in  Figure  1(e)  we  show  the  computed  corresponding  region,  shown  in  red  on 
the  outer  hull  surface.  Figure  2  illustrates  the  axial,  coronal,  and  sagittal  slices  of  the  original 
MRl  brain  volume  on  which  we  overlay  the  intersection  curve  of  the  pial  and  outer  hull  surfaces 
of  one  hemisphere  (in  blue)  with  the  selected  sulcal  region  on  the  pial  surface  and  the  computed 
corresponding  region  on  the  outer  hull  surface  (in  red). 

These  geometric  computations  form  the  basis  needed  for  the  gyrification  indexes 
presented  next.  Before  that,  we  should  note  that  in  contrast  with  (Schaerei  al.  2008),  where  the 
intersection  of  an  extrinsic  sphere  of  predefined  radius  with  the  outer  hall  is  used,  our 
computation  of  the  corresponding  outer  hull  region  is  completely  intrinsic  to  the  geometry  of  the 
prc^lem. 
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2.3  (ryrification  Indexes 


The  most  commonly  used  gyrification  index  was  proposed  by  Zilles  et.  al.  in  1988  (ZiUes 
el  al.  1988)>  and  is  defined  as  (he  ratio  between  the  lengths  of  the  total  cortical  surface  and  the 
superficially  exposed  cortical  surface  on  a  two  dimensional  slice  in  the  coronal  section, 

where  is  the  length  of  the  total  cortical  surface  on  a  slice  > .  and  is  the  length  of  the 
superficially  exposed  cortical  surface  on  the  slice  t .  The  mean  GI  of  each  hemisphere  is  then 
defined  as 


01 


where  the  sum  runs  over  all  the  slices  t .  This  measuremeni  does  no(  allow  for  the  assessment  of 
localized  properties,  including  development,  of  the  gyrification.  Different  theories  have 
postulated  that  active,  localized  growth  of  the  cortical  convolutions  may  lead  (o  different  brain 
siruciures.  e.g.,  regions  which  are  developing  into  gyri  grow  at  a  faster  rate  than  areas  destined  (o 
be  become  sulci.  Furthermore,  recent  research  demonstrates  ihai  the  development  of  differeni 
regions  of  the  brain  may  result  in  remarkably  different  brain  functioning,  e.g,,  the  temporal  lobe 
is  more  vulnerable  to  insults  from  preterm  births,  with  a  resulting  increased  temporal  lobe  GI 
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(KcAler  ei  al.  2006).  Thu^.  reliable  and  accurate  localized  gyrification  indexes  are  of  great 
intere.s(  in  the  understanding  of  (he  complexity  of  brain  folding. 


A.  Arta-Bostd  Gyrificalion  Meaxuremenl 

We  now  present  the  first  proposed  local  gyrificalion  index  measuremenL  ,  Using 
freesurfer  (Dale  etal.  1999).  each  hemisphere  is  classified  into  35  different  sulcal  and  gyral 
regions,  and  a  complete  labeling  of  these  35  regions  can  be  obtained  by  the  auiomated 
parcellaiion  system.  We  group  these  ROIs  into  six  regions  (see  Table  1)  of  (he  human  cortex: 
fronial  lobe,  parietal  lobe>  temporal  lobe,  medial  lemporal  lobe,  occipital  lobe,  and  the  cingulaie 
gyrus. 

Figure  3  shows  the  classification  of  the  6  regions  on  (he  left  pial  surface.  We  can 
successfully  identify  the  corresponding  6  regions  on  the  hull  surface.  For  comparison >  we  also 
show  the  inflated  left  pial  surface,  which  we  obtain  directly  from  freesurfer  (Dale  et  al.  1999). 

Our  local  gyrificalion  index  Cil  ^  at  each  lobe  is  (hen  defined  as  the  ratio  between  the 
area  of  the  lobe  on  the  pial  surface  and  the  area  of  the  corresponding  region  on  the  outer  hull 


surface. 
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where  A’  is  the  area  of  any  selected  lobe  on  ihe  pial  surface  and  a'  is  the  corresponding  area  on 
ihe  hull  surface,  computed  as  detailed  in  the  previous  section. 

This  local  gyrification  index  (U  ^  has  several  key  features: 

I .  The  corresponding  outer  hull  area  is  naturally  computed  from  the  ROI  on  the  pial  surface, 
leading  to  a  parameter  free  intrinsic  GI,  If  we  select  the  whole  pial  surface,  the 
corresponding  region  is  then  the  whole  outer  hull  surface,  and  the  gyrification  index 
measurement  becomes 


where  A ,  is  the  area  of  the  whole  pial  surface,  and  a.  is  the  area  of  the  whole  outer  hull 
surface.  If  we  only  consider  a  slice  of  the  pial  surface,  the  surface  area  in  3D  becomes  the 
slice  perimeter  in  2D,  However,  this  reduction  may  still  be  different  from  the  exact 
definition  of  Zilles  (Zilles  et  al,  198^).  since  the  corresponding  point  on  the  outer  hull 
surface  may  not  be  on  the  same  plane  any  longer.  From  this  point  of  view,  powerful  3D 
gyrification  index  computation  is  of  great  importance. 

2.  Since  all  the  triangles  in  the  selected  region  on  the  pial  surface  are  collected,  all  the  sulcal 
information  restricted  to  the  selected  area  is  considered. 

3.  The  selection  of  the  region  on  the  pial  surface  can  be  further  parceled  to  determine  Gl  for 
different  subregions.  For  instance,  having  calculated  the  depth  information  for  each  point 
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on  ihe  pial  surface,  we  can  define  a  threshold  and  choose  a  region  of  interest  as  aJl  ibe 
points  with  depth  rf  a  4/^. 


b.  Dcpih'BMed  Cvrf/icc/tft»n  MeusuTemtnt 

Before  we  extend  the  gyrification  index  QI  to  a  weighted  gyrification  index  by  adding 
the  depth  information  of  the  piaJ  surface,  we  introduce  a  normalization  of  the  computed  depth 
value. 

Since  the  absolute  depth  significantly  depends  on  brain  size,  it  is  reasonable  to  consider 
the  relative  depth.  In  this  way,  two  brains  with  a  similar  morphology,  but  different  volumes, 
would  have  comparable  gyrification  index.  We  implement  this  by  using  rf,y  instead  of  the 
absolute  value  d . 

j 

d.  ■ - . 

yy  !A 

where  the  denominator  iv  /  a  on  the  right  hand  side  is  a  normalization  factor  based  on  volume 
v  and  area  A  {Rodriguez-Carranza  et  aJ.  2008).  We  can  examine  this  normalization  considering 
a  sphere  of  radius  R  ,  then  the  normalized  value  for  a  certain  depth  d  of  folding  becomes 

rf  J 

d,  s  — - s  — . 

3 


and  if  the  radius  is  and  the  depth  is  .  the  normalized  depth  becomes 
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2J  d 

J,  »  — - «  - 

3 

While  these  two  spheres  have  similar  folding  but  different  depths,  the  normalized  depths  are  the 
same. 

We  now  define  the  depth-based  gyrification  index.  For  any  region  selected  on  the  piaJ 
surface,  the  corresponding  region  on  the  outer  hull  surface  is  denoted  by^.,  A  gyrification  index 
that  considers  the  normalized  depth  information  can  be  defined  as 


G( 


k>c^ 


where  and  M,  are  the  total  number  of  triangular  vertices  inside  the  regions  and  R,, 
respectively,  d[.  is  the  normalized  depth  of  the  vertex  t  (depth  to  the  outer  hull  surface  as 
detailed  before  (Kao  et  aJ.  2007),  is  the  area  of  the  vertex  •  in  ,  and  4/  is  the  area  of  the 
vertex  j  in  .  Our  definition  of  the  area  at  a  vertex  •  (see  Figure  4.  cyan  colored  vertex)  is 
defined  in  the  following  way.  For  each  neighboring  triangle  of  vertex  t ,  we  find  the  centroid  (see 
Figure  4,  black  colored)  and  for  each  neighboring  ring  of  vertex  r ,  find  the  middle  point  (see 
Figure  4,  yellow  colored).  Then  the  area  at  vertex  i  is  the  area  sum  of  all  such  triangles  formed 
by  the  vertex  t .  the  centroid  of  the  neighboring  triangle  and  the  middle  point  of  the  neighboring 
ring  in  the  same  triangle.  Figure  4  illustrates  this. 
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By  adding  weights  lo  (he  depth  of  ihe  pial  surface  on  (he  numerator  of  Gi  ^  .  we  are 
approKimadng  the  volume  of  ihe  res(ric(ed  par(  be(ween  (he  selected  region  on  ihe  pial  surface 
and  the  corresponding  region  on  (he  ou(er  hull  surface.  Thus,  the  3D  local  gyrification  Index 
Cl  is  in  fact  a  measurement  of  the  average  geodesic  depth  in  the  selected  region  of  Interest. 
Similar  lo  the  previously  defined  gyrification  index,  the  global  counierpan  of  GI  ^  becomes 


Of 


i 


S''.’ 


This  is  an  approximation  of  the  average  geodesic  depth  of  the  restricted  (CSF)  region  between 
the  outer  hull  and  the  pial  surface  { Y  s  o  and  cp  2  Oj , 


C.  Pointwi.^e  Gyrification  Mtasuremenl 

Based  on  the  geodesic  depth.  Section  2.1,  we  now  present  a  point-wise  GI.  That  is,  for 
any  vertex  >  on  (he  pial  sur^e.  with  its  corresponding  normalized  geode.sic  depth  rf' .  we 


define 
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There  is  a  connection  between  this  definition  and  ,  let  us  see  this  connection  now. 
Let  us  assume  that  the  selected  region  of  interest  on  the  pial  surface  has  total  area  <4^ .  and  the 
corresponding  region  on  the  outer  hull  surface  has  total  area  4/ ,  Then. 


<31 L 


«»i _ 

If. 

Z''.' 

>»i _ 

z 


*Z‘';-^;'Z'^; 


<•1 


/■I  lO 

_ f 

A' 


-  A' 


m  d. 


where  d ,  is  the  average  geodesic  depth  in  the  selected  region  on  the  pial  surface.  As  the  selected 
region  gets  smaller  and  smaller,  and  finally  converges  to  the  vertex  t ,  it  becomes  natural  to 
define  the  local  gyrification  index  at  the  vertex  using  its  own  normalized  geodesic  depth, 
obtaining  GI  ^  .  {Note  that  we  dropped  the  ratio  between  4 '  and  4.'  foe  simplicity ) 
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2.4  ComparlsAD  with  previous  work 


In  comparison  with  (Schaer  ei  al.  2008),  our  proposed  G1  compuiaiion  method  1% 
parameter  free  and  defines  the  local  3D  gyrification  index  directly  on  the  pial  surface.  In  {Schaer 
el  al.  2008),  they  redi.siribuie  the  local  gyrification  index  to  any  vertex  point  v,  on  the  pial 
surface  by  first  defining  the  local  gyrification  index  on  the  outer  hull  surface,  and  the  weighting 
during  redistributing  Is  related  to  the  distance  between  v,  and  the  involved  vertex  points  ,  f  on 
the  outer  hull  surface.  There  are  several  possible  disadvantages  of  this  method:  first  of  all.  the 
outer  hull  surface  is  a  surface  generated  from  the  pial  surface,  and  .strictly  speaking,  there  U  no 
gyrification  index  on  the  outer  hull  surface.  Secondly,  the  extrinsic  distance  they  use  in  the 
weighting  is  along  the  normal  axis  of  the  vertex  to  the  outer  hull  surface,  while  our  intrinsic 
distance  is  the  computed  geodesic  depth,  which  is  a  more  natural  technique  in  defining  the 
gyrification  index  (Kaoet  al,  2007),  In  addition,  after  Schaer  et  al.  (Schaer  et  al.  2008)  extract 
the  perimeter  on  the  region  of  interest  on  the  outer  hull  surface,  they  project  this  perimeter  to  the 
pial  surface  by  taking  the  nearest  point  on  the  pial  surface,  whereas  our  correspondence  relies 
once  again  on  the  natural  geodesic  depth  direction.  Finally,  the  results  from  a  spherical  mask  will 
depend  on  the  radius  used.  Smaller  spheres  will  have  greater  variability  secondary  to  anatomical 
differences,  whereas  larger  spheres  will  not  provide  a  localized  enough  measure  of  gyrification 
and  will  cross  boundaries  into  alternate  regions.  Nevertheless,  the  approach  by  Schaer  et  al. 
(2008)  U  a  significant  improvement  over  the  previous  2-0  Gl. 


Geometric  Corrputation  of  Gyrification  Indeaes 


ChiU'Ven  Kao 


20 


2.5  Subjects 


To  study  ihe  development  of  gyrification,  we  implement  our  method  on  a  sample  of  26 
typically  developing  children  and  adolescents.  Subjects  were  stratified  into  three  age  groups, 
which  included  children  to  12  years  old),  young  adolescents  (12  to  15  years  old),  and  older 

adolescents  ( 1 6  to  19  years  of  age).  The  demographic  information  for  the  children  is  provided  in 
Table  2.  Socioeconomic  status  was  determined  using  the  Four  Factor  Index  of  Social  Status 
(Hollingshead,  1975). 


2.6  Statistical  Analyses: 


Evaluation  of  the  differences  in  demographics  wa.s  assessed  using  either  one* way 
ANOVA  or  t-tests  for  continuous  data  and  chi-square  le.sts  for  categorical  data.  Paired  t-tesLs 
were  used  to  assess  lateralization  of  gyrification  and  cortical  depth  measures.  Mea.sures  which 
did  not  demonstrate  lateralizing  effects  were  pooled.  An  ANCOVA  of  the  total  sample  (14  males 
and  12  females),  with  age  and  ICV  as  acovariate,  was  used  to  evaluate  sex  effects  of  sulcal 
depth.  Measures  of  gyrification  were  evaluated  between  the  three  age  groups  using  a  one-way 


ANOVA. 
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3  Exp«rimen(al  R^ults 


3.1  Demographics 


There  were  no  differences  between  the  sex.  socioeconomic  status,  or  handedness  between 
any  of  the  three  age  groups  (T^le  2),  In  addition,  when  stratified  into  two  groups  by  sex.  there 
were  no  differences  in  the  demographic  data.  An  ANCOVA  with  age  as  a  covariaie 
demonstrated  sex  differences  in  the  total  intracranial  volume  (Fi  jj »  4,8,  p  s  0.04).  however, 
there  were  no  differences  in  ICV  between  the  three  age  groups. 


3.2  Data  Aoaly'sis  on  the  Pial  Surface 


Before  analyzing  the  data  in  detail,  we  present  an  overview  of  the  area  ratio  and  average 
normalized  depth  for  each  Jobe  on  the  pial  surface.  The  six  regions  are  shown  in  Figure  3, 
marked  with  different  colors.  For  simplicity,  in  the  following  tables,  we  denote  the  frontal  lobe 
by  'F*.  the  parietal  lobe  by  the  temporal  lobe  by  ’T',  the  medial  temporal  lobe  by  MT',  the 
occipital  lobe  by  ’O’,  and  the  cingulate  cortex  by  'C’. 
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Figure  5(a)  plots  the  area  ratio  —  and  Figure  5(b)  plots  the  normalized  average  depth 

A 

I* 


for  each  lobe  in  l>oih  hemispheres  for  subjects  in  the  control  group.  Here,  ibe  notations 


are  the  same  as  we  mentioned  in  Section  2  and  the  labeling  names  and  the  grouping  of  the  sulci 
for  each  lobe  are  listed  in  Table  1 .  Tables  2  and  4  list  the  mean  value  of  the  area  ratio  and  the 
normalized  average  depth  for  the  control  group.  They  both  show  no  significant  inirin.sic 
differences  in  the  left  and  right  hemispheres. 


3.3  Data  Analysis  of  the  Outer  Hull  Surface 


The  area  ratio  —  for  the  corresponding  region  of  each  lobe  on  the  outer  hull  surface  are 

shown  in  Figure  6  and  Table  5.  As  expected,  there  are  no  significant  differences  between  the 
hemispheres  and  the  corresponding  ordering  is  as  for  the  pial  surface.  This  supports  the 
technique  we  designed  to  find  correspondences  between  the  pial  surface  and  the  outer  hull. 
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3.4  Lateraliialion  of  Sulcal  Depth  aod  (ryrification  Measures 

A  paired  t-(esi  found  (hat  (here  was  significant  laterality  for  the  maximum  depth  for  many 
of  the  major  sulci.  These  iaierabty  difference.^  are  shown  in  Table  6.  As  expected,  the  regions 
which  demonstrate  iateralizing  effects  are  related  to  language  regions  in  the  temporal  lobes,  as 
well  a.s  sulci  within  the  frontal  and  parietal  lobes.  Sulci  that  do  not  show  Iateralizing  effecLs 
include  areas  related  to  sensory  and  motor  functions  (central  sulcus  and  calcarine  fissure).  Due  to 
these  differences,  lateralized  measures  were  not  pooled  to  assess  for  developmental  differences 
in  sulcal  depth. 

Paired  t-tests  found  no  differences  in  either  Gl  ^  or  GI  ^  between  the  left  and  right 
frontal  and  temporal  lobes.  However,  there  were  significant  laterality  effects  of  both  Gl  ^  and 

in  the  parietal  ld3e(Gi;^^  :  t  s  4.8.  df*  28.  p  <  0.0001;  GI^  :  t  s  2.9.  df  s  28.  p  s 
0.007),  occipital  lobe  (Gl  :  t  s  4.6.  df  s  28.  p  <0.0001;  :  t  s  4.6,  df  s  28,  p  s  0.01).  and 

the  cingulate  cortex  (Gl|^  :  t»3,3,dfa28.  p  a  0,004;  CI^  :la  2,2,  df  a  28,  pa  0.03).  The 

parietal  lobe  had  greater  cortical  complexity  on  the  right,  whereas  the  occipital  lobe  and 
cingulate  cortex  had  greater  cortical  complexity  on  the  left.  In  addition,  laterabiy  effects  were 
found  in  the  medial  temporal  Id^es  for  d  ^  (t » *2.4,  df  a  28,  p  a  0.02).  with  the  left  lobe 
having  greater  gyrification  than  the  right.  Due  to  these  differences,  gyrification  measures  were 
not  pooled  for  the  right  and  left  hemispheres. 
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3.5  Devdopmental  Differences  in  Sulcal  Depth 

There  were  reJaiively  few  differences  in  sulcal  depth  between  the  different  age  groups  for 
the  15  major  sulci  that  were  evaluated.  There  were  age-related  differences  in  the  maximum  depth 
of  the  right  middle  temporal  suicua  (P221  s  4, 1 ,  p  s  0.03).  right  lateral  sulcus  (F221 »  3.5.  p  » 
0.05).  and  a  trend  for  the  right  central  sulcus  (Fj  jj  *  3.2.  p  s  0.06), 


3.6  Developmental  Differences  in  (>yrification 


There  were  no  differences  in  the  mea.sures  of  Gi  or  Gi  ^  between  males  and 
females  in  any  of  the  regions  of  intere.st.  thus  these  measures  were  pooled  for  the  analyse,s.  Age- 
related  differences  in  gyrification  for  01  in  the  left  and  right  hemispheres  can  be  seen  in 
Figure  7(a)  and  Figure  7(b),  respectively.  The  GI  decreases  with  age  in  the  left  (Fj  jj*  4,6,  p 
s  0,02)  and  right  (F1.24S  4.9.  p  s  0.02)  frontal  lobe;  left  IF224  *  4.2.  p  »  0.03)  and  right  (F;24® 
6.8,  p  s  0.004)  medial  temporal  lobe.  In  addition,  there  are  trend  decreases  in  the  GI  1^,  by  age 
in  the  left  (Fi.m  s  2.6.  p  »  0.09)  and  right  parietal  (F22a »  2.6,  p  s  0.09). 

Fewer  age-related  differences  were  found  u.sing  the  gyrification  measure  which  accounts 
for  regional  differences  in  sulcal  depth  (GI  ^  ).  Age-related  differences  were  found  in  the  left 
(F’.id®  4.1.  p  s  0.03)  and  right  (Fj  jj®  6,5,  p  s  0.006)  medial  temporal  lobe.  In  addition,  trend 
differences  were  found  in  the  left  (F’  ld  s  3. 1.  p  s  0.07)  and  right  (Fj.y  »  2,8,  p  s  0.08)  cingulate. 
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left  parieiaJ  lobe  (F2  24»  2,7,  p  s  0.09);  and  ihe  left  (F2.’4  s  2.5.  p  »  0.10)  and  right  (F2,24s  2.3.  p 
s  0,06)  frontal 


3.7  Comparison  with  Existing  Techniques 


To  compare  our  technique  with  ihe  local  gyrification  index  (IGI)  proposed  by  Schaerel  al. 
(Schaer  et  ai,  200^),  we  calculated  iGl  results  at  different  radii  and  averaged  these  results  into  the 
six  defined  brain  regions.  It  turns  out  that  IGI  is  quite  sensitive  to  the  selection  of  the  radius  of 
the  sphere,  {Radius  of  20, 25,  30.  and  35nini  as  proposed  are  tested.)  With  radii  of  less  than 
20mm.  not  all  data  could  be  processed  using  Frees urfer  (Dale  el  al.  1999)  with  a  resultant  los.s  in 
the  degrees  of  freedom.  Table  7  lists  the  one*way  ANOVA  results  for  radius  of  20,  25,  30.  and 
35mm  respectively. 

To  illustrate  this  strong  radii  dependency,  we  ran  a  one-way  ANOVA  test  (Table  8)  to 
show  the  sen.sitivity.  From  Table  8.  we  can  observe  significant  differences  between  radii  in  the 
frontal,  parietal,  and  temporal  lobes  in  both  hemispheres  for  the  three  age  groups.  More 
variability  appears  if  we  do  analyses  using  gender  as  a  covariate  while  the  measures  Gi  and 
Gi  ^  we  proposed  behave  more  stable.  To  have  a  more  complete  idea  of  the  effects  of  IGI.  we 

did  a  four-radius  by  three-age-group  by  rwo«region  ANOVA  test,  which  showed  that  there  were 
significant  effects  in  radius  (F3,i27i  s  44.61.  p  <  0.0001).  age  group  (F2.i27i  s  19.03.  p  <  0.0001) 
and  no  significance  in  the  region  side  (Fm27)  s  l.I  I.  p  s  0,29). 
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3.8  Discussioo 


We  utilized  &  novel  gyrification  index  computational  technique  to  demonstrate 
developmental  differences  in  gyrification  in  a  cohort  of  typically  developing  children  and 
adolescents.  The  proposed  computational  method  employed  intrinsic  3D  measurement  to  find 
gyrification  indexes  for  a  region  of  interest  or  a  single  point  on  the  pial  surface.  It  i.s  thus  able  to 
identify  local  alternations  in  brain  surface  complexity. 

Our  results  support  continued  molding  of  the  surface  morphology  of  the  brain  through 
childhood  and  adolescence.  This  same  time  period  has  been  shown  to  be  associated  with  a 
decrease  in  surface  GM  and  an  associated  increase  in  the  surface  cerebral  spinal  fluid  (Giedd  et 
al.  1999;  Sowell  et  al.  2004).  It  is  not  surprising  that  decreases  in  the  underlying  volume  of  GM, 
would  result  in  an  associated  decrease  in  conical  complexity,  as  has  been  demonstrated  in 
typically  developing  adults  (Magnotta  et  al.  1999).  The  developmental  differences  in  GI  are 
localized  to  the  later  maturing  structures,  including  the  frontal  and  medial  temporal  Irt^es 
(Huttenlocher,  1982. 1990).  Regions  that  develop  earlier,  such  as  the  occipital,  parietal,  and 
lateral  temporal  lobe.s,  do  not  demonstrate  the  same  level  of  developmental  differences  as  the 
frontal  and  medial  temporal  regions  (Huttenlocher  and  Courten.  19S7). 

The  brain  has  reached  its  peak  volume  by  eight  years  of  age  and  thus  it  is  not  surprising 
that  there  are  not  striking  developmental  differences  in  the  depths  of  the  major  sulci.  Indeed,  the 
primary  sulci  are  the  first  to  develop  and  are  associated  with  a  greater  herit^ility  (Lohmann  et  al. 
1999).  The  secondary  and  tertiary  sulci  develop  later,  and  are  more  linked  to  environmental 
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factors.  However,  most  of  the  primary  and  secondary  sulci  form  during  the  third  trimester  of 
fetal  life  (Welker.  1990)  and  it  is  likely  that  restructuring  and  depth  of  the  major  sulci,  including 
laterality  differences,  are  fairly  stable  by  eight  years  of  age.  Future  studies  mapping  the 
development  of  the  sulci,  gyri  and  G1  between  birth  and  eight  years  of  age  will  be  beneficial  to 
assess  alterations  of  the.se  measures  as  language  and  other  cognitive  processes  mature. 

There  are  limitations  to  the  current  study.  We  had  a  relatively  small  sample  size  to 
evaluate  for  developmental  differences.  However,  the  fact  that  significant  findings  were  present 
and  were  consistent  across  measures  may  reflect  a  relatively  robust  measure  to  assess 
developmental  differences. 


4  Concluding  Remarks  and  Discussion 


Our  3D  geometric  approach  for  defining  the  global  and  local  gyrification  indexes 
naturally  extends  the  standard  way  of  defining  the  coronal  2D  GI  and  fully  utilizes  the  intrinsic 
3D  nature  of  human  cortical  surface.  Classifying  the  cortical  surface  into  6  different  regions  and 
defining  a  simple  local  gyrification  index  on  each  region  as  the  area  ratio  of  corresponding 
surfaces,  helps  us  to  more  carefully  characterize  the  cortical  complexity  in  each  region  and 
therefore  provides  a  better  way  to  understand  the  effect  of  each  region  in  brain  functioning. 
Furthermore,  applying  our  method  in  a  clinical  study,  the  proposed  measurement  of  depth- 
weighted  local  gyrification  turns  out  to  be  more  robust  in  finding  the  developmental  differences 
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in  children  and  adolei^cenls.  e.g,>  we  observe  significaniJy  increased  gyrification  in  the  right 
parietal  lobe  and  right  cingulate  cortex,  as  well  as  age*rated  differences  in  the  left  frontal,  right 
parietal  and  the  right  cingulate  cortex.  These  findings  provide  some  references  for  future  study 
of  the  relationship  between  gyrification  and  neurological  and  psychiatric  conditions,  in  addition 
to  the  development  of  other  more  advanced  techniques  to  quantify  the  gyrification  of  the  human 
brain. 

Since  there  is  not  one  universally  accepted  definition  for  computational  GI,  and  this  is  a 
relatively  new  and  ongoing  research  area,  we  conclude  this  discussion  with  some  possible 
alternatives  that  deserve  further  study. 

From  a  statistical  point  of  view,  the  3D  local  gyrification  index  GI  ^  is  essentially  an 
estimation  of  the  first  moment  of  the  area.  This  motivates  us  to  consider  higher  order  moments. 
For  in.stance.  a  second  central  moment  estimator 


GI 


* 


where 


do 


<  I 
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is  the  average  depth  in  the  selected  region  and>  evaluates  ibe  variance  of  ihe  area  distribution  and 
the  third  central  moment  estimator  will  then  help  in  evaluating  the  skewness  of  the  area 
distribution. 

The  gyrification  index  measurement  does  not  necessarily  need  to  be  scalar  valued,  we  can 
also  consider  a  vectorized  gyrification  index  measurement.  For  any  selected  region  on  the 
pial  surface  and  the  corresponding  region  on  the  outer  hull  surface,  we  define  (for  example) 


GI 


MoJ 


GI 


fefW 


fGlL^.GlL  I 


A 

_ f 

A. 


1^.-' 


H. 


With  this  expression,  we  can  evaluate  the  local  sulci  information  considering  both  the  area  ratio 
and  the  average  depth  of  the  selected  region,  narrowing  the  similarity  between  different  regions 
and  providing  more  information  for  analyzing  the  brain  complexity. 

These  alternatives  for  GI  computations  will  be  tested  and  reported  elsewhere. 
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Table  1:  Parcellation  of  (he  35  regions  into  6  regions. 


lobe  name 

region  names 

Formal 

caudal  middlelronial.  laieraJorbiioln^ntal.  medialorhiiofronial. 

paracentral.  parsopercularis>  parsorbitaJis.  parstnangularis. 

precentral.  rosiraJmiddlefrontal>  superiorfromal,  fronialpole 

Parietal 

inferiorparietal.  p<»tceniraJ.  precuneus,  superiorparieial. 

supramarginal 

Temporal 

bankssts.  inferiortemporaJ.  middletemporal.  siiperionemporal. 

temporalpole.  (ransversetemporal 

Medial 

Temporal 

enloriiinal.  fusiform,  parahippocampal 

OccipitaJ 

cuneus,  lateraloccipitai,  lingual.  pericaJcarine 

Cingulate 

caudal  anteriorcingu  late,  i^thm  use  ingu  late,  posieriorcingulate. 

rosiraJanleriorc  i  ngulate 
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Table  2:  Demo^aphic  Information  for  the  26  Typically  Developing  Children  between  (he 
ages  of  S  and  19  Years  of  Age. 


8  to  12  Year 
Olds 

13  to  15  Year 
Olds 

16  to  19  Year 
Olds 

P 

Age 

10.5  (1.6) 

i4U  tU8?) 

16.9  (1.2) 

Sex  (Male  / 
Female) 

4/2 

6/3 

4/7 

n.s. 

Socioeconomic 

Status 

50.7  (9.7) 

49.7(15.5) 

51.1  (9.7) 

n.s. 

Table  3:  Average  area  ratio  of  each  region  for  subjects  in  (he  control  group.  'Ih'  stands  for 
(he  left  hemisphere  and  'rh'  stands  for  rbe  right  hemisphere. 


Kcbicm 

i- 

K 

1 

M*I* 

() 

C 

Ih 

0.3509 

0,2265 

0.1540 

0.0484 

0,1125 

0.0364 

rh 

0.3533 

0,2302 

0.1482 

0.0441 

0,1164 

0,0380 
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Table  4:  Meao  value  of  the  DormaJized  average  depth  of  each  region  for  subjects  in  the 
control  group.  'Ih'  stands  for  the  left  hemisphere  and  'rh*  stands  for  the  right  hemisphere. 


Kcvicm 

i- 

P 

1 

M*I* 

() 

C 

Ih 

0.1241 

0,1613 

0.1501 

0.0800 

0,1001 

0,0637 

rh 

0.1230 

0,1625 

0.1507 

0.0757 

0,0986 

0.0649 

Table  Si  Average  area  ratio  of  each  corresponding  region  of  each  hemisphere  for  subjects 
in  the  control  group.  'Ih'  stands  for  the  left  hemisphere  and  'rh'  stands  for  the  right 
hemisphere. 


Kcvidii 

i' 

P 

1 

M*I* 

() 

C 

Ih 

0.3750 

0  1937 

0  1527 

0.0634 

0.1180 

0  0509 

rh 

0.3777 

0,1909 

0.14S9 

0.0598 

0,1261 

0,0507 
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Table  6:  Maximum  depth  (in  mm)  for  the  major  suJci  by  hemUphere  and  paired  ('test 
showing  lateralization. 


Suki 

Maximum 
Depth  (left) 

Maximum 

Depth 

(right) 

t/p 

C'oiitrjl  Sukus 

30.1  (2.4) 

n.s. 

Ljlcral  Sulcus 

.14  1  (2  1) 

33.1(1  9) 

2.5/0  02 

Supcriitr  iTunlal 

20  M  1  6) 

192  (2  2) 

2.1/0.05 

Superior  1cm  pordi 

22f.(22) 

24.8  (2.2) 

-7  5/ 
<0,0001 

Middle  Temporal 

14.2(1.9) 

1S.5(2.I) 

-2.6/0.02 

Inlrupaiielal 

260(1  7) 

26.0(2.1) 

n.s. 

Transverse  Occipital 

172(1.9) 

16.4  (2.3) 

n.s. 

liliiHH 

20  2(3  2) 

20.8  (3.5) 

n.s. 

Superior  Parietal 

13.3(1.7) 

14.3(1.3) 

-2.9/0  008 

Pariet«>occipiial 

Fissure 

18  8(2.6) 

20.5  (2.6) 

Supramarginal 

21  9(22) 

21.1  (3.6) 

n.s. 

Orbilofrontal 

9.6  (1.3) 

10.8(1.0) 

-4  7/ 
<0,0001 

Ant.  Ascending 

Ramus 

29.7(2.8) 

27.3  (2.2) 

4.8/ <0.0001 

Ant.  HorizoniaJ 

Ramus 

13.3(1.7) 

16.2(3.3) 

18.5/ 

<0,0001 
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Table  7:  One-way  ANOVA  test  show  log  develApmeDta)  dilTerences  between  IGI  for 
different  radius  Imm)  within  different  regions. 


Brain 

Region 

LGI  rs2A 

L(;l  rs25 

L(;l  rs30 

L(;l  rs35 

Left 

ftjVp 

■^21 

Uft 

Fiu/p 

Right 

Fj24/P 

Left 

F2ja/p 

Right 

FjWp 

Uft 

Fiu/p 

Right 

F224/P 

i'O'nlj) 

5  ?/0()l 

.V I/O  07 

2  7/0  (N 

2  5/(1  10 

4A).05 

2.VA)07 

7A)  04 

^  0/0  07 

Panoldl 

I/O  (Mil 

5.2A)OI 

7  7/0  00'^ 

0/(1  07 

t  I/O  (Ml7 

4.4A)02 

6  .<A)  (H)6 

4A)()^ 

1  ciiiporjl 

5  I/O  01 

4.5  A)  02 

2/0  (lf> 

2  5/0  10 

2  6A).10 

.V4A)05 

1  9A).2 

^  .Vo  06 

Mediul 

Temporal 

0  4/0  7 

4.5/0  02 

0  .m? 

m 

0.5A)  h 

1  5/0  2 

0  .mx 

I.IA).^ 

Occ  ipiul 

2  2/0  1 

1  5/(1 

^am 

■Kia 

2  Vrti.o? 

1  2/0 

2  7A)  (W 

2.1A)  1 

Cioyuldlc 

5  X/0(MW 

.V7/0  04 

4  7/0  (12 

4  fi/0  02 

4 

5.6A)OI 

4  M)  02 

4  .Vo  0.^ 

Table  8:  One-way  ANOVA  test  showing  differences  of  IGI  between  radius  (mm)  for  each 
age  group  in  each  region. 


Brain 

Region 

Agegroupif^l2i 

Agegrnupi  13-15 1 

Agegrnupi  16-19) 

■SSB 

I'nmUil 

9  1/0.0005 

1 1  5/0.0001 

7  7A)0004 

>i.lAKK)03 

6  .5/0  (Mil 

10  9/0 
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Figure  1:  Flow  chart  of  (he  main  steps  of  our  proposed  GI  compulation  algorithm.  Each 
figure  explains  the  major  steps  of  our  algorithm: 


(a)  Pial  surl'ace  extraction 


(bi  Outer  hull  surface  computation 
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(e)  Geodesic  depth  compuLniion 


(dj  Region  of  interest  on  the  pial  surface  selection 
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(e)  Desponding  region  on  ihe  outer  hull  surface  compulation 
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Figure  2:  AxJaJ.  coronal,  and  fmgittal  slices  of  the  MRI  brain  volume  combined  with  the 
(dul  and  outer  hull  surfaces  of  one  hemis^there  (blue  lines).  The  red  pari  on  the  curve  is  the 
intersection  of  the  selected  region  {frontal  lobe  here)  on  the  pial  surface  und  the 
corresponding  region  on  the  outer  hull  surface.  Top  three  slices  correspond  to  Figure  1(d) 
(pial  surface)  and  bottom  three  slices  correspond  to  Figure  1(e)  {outer  hull). 


Geometric  Computation  of  Gyrification  Indexes 


ChiU'Ven  Kao 


42 


Figure  3:  6  on  left  pial  surface,  InHaicd  left  pial  surface,  and  the  left  hull  surface. 

Frontal  lube  is  shown  in  blue,  parietal  lobe  in  green,  medio  temporal  lobe  in  cyan,  occipital 
lobe  in  red.  and  cingulaied  cortex  In  magenta.  The  region  which  is  not  classified  is  shown  in 
brown. 
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Figure  4:  Area  derinition  a(  a  vertex.  For  the  vertex  at  the  center  (c^an).  the  area  Is  defined 
as  the  ivtal  area  of  the  region  !>hown  in  magenta,  which  Is  made  up  uf  aJI  the  triangles 
constructed  by  connecting  the  vertex  (cyan),  the  centroid  I  black  I  of  its  neighboring  triangle 
and  the  middle  point  (yellow)  of  its  neighboring  ring  in  the  same  triangle. 
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Figure  5(a):  Area  ratio  of  each  region  on  both  hemispheres  for  subjects  in  the  control 
group.  Left  beioisphere  measuremenu  are  marked  by  and  right  hemisphere 
measureroents  are  marked  bv 


Figure  5(bi:  Normalized  average  depth  of  each  region  on  both  hemisplieres  for  all  subjects. 
Left  hemisphere  measurements  are  marked  by  *+*.  and  right  hemisphere  measurements 
are  marked  by  *•'. 
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Figure  6:  Area  ratio  of  the  corresponding  region  for  each  region  on  both  hemispheres  for 
ail  subjects.  Left  bemispbere  measurements  are  marked  by  and  right  beiuispherv 
measureroents  are  marked  bv 
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Figure  7:  DvveJopmujiiai  differoDccs  bet^veeo  gyrification  indices  within  different  regions 
fur  the  both  hemispheres. 
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